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ABSTRACT 

This is the second in a series of papers presenting results from the SCUBA Local Uni- 
verse Galaxy Survey. In our first paper we provided 850/iin flux densities for 104 galax- 
ies selected from the IRAS Bright Galaxy Sample and we found that the 60, 100/im 
(IRAS) and 850^m (SCUBA) fluxes could be adequately fitted by emission from dust 
at a single temperature. In this paper we present 450/im data for the galaxies. With 
the new data, the spectral energy distributions of the galaxies can no longer be fitted 
with an isothermal dust model - two temperature components are now required. Using 
our 450/xm data and fluxes from the literature, we flnd that the 450/im/850/xm flux 
ratio for the galaxies is remarkably constant and this holds from objects in which the 
star formation rate is similar to our own Galaxy, to ultraluminous infrared galaxies 
(ULIRGS) such as Arp 220. The only possible explanation for this is if the dust emis- 
sivity index for all of the galaxies is ~ 2 and the cold dust component has a similar 
temperature in all galaxies [Tc ~ 20 — 21 K). The 60/im luminosities of the galaxies 
were found to depend on both the dust mass and the relative amount of energy in 
the warm component, with a tendency for the temperature effects to dominate at the 
highest Lgo- The dust masses estimated using the new temperatures are higher by a 
factor ~ 2 than those determined previously using a single temperature. This brings 
the gas-to-dust ratios of the IRAS galaxies into agreement with those of the Milky 
Way and other spiral galaxies which have been intensively studied in the submm. 
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1 INTRODUCTION 

This paper will present further results from the SCUBA Lo- 
cal Universe Galaxy Survey (SLUGS), which is the first sys- 
tematic submillimetre survey of the local universe. The gen- 
eral aims of the survey are to provide the statistical measure- 
ments (the submillimetre luminosity and dust-mass func- 
tions) necessary for interpreting observations of the high- 
redshift universe and to determine the submm properties of 
dust in the local universe. The ideal way to carry out the lo- 
cal survey would be to do a blank field survey of large parts 
of the sky but this is currently impractical because the field 
of view of SCUBA is only ~2 arcmins. Our approach was 
to observe galaxies drawn from as many different complete 
samples selected in as many different wavebands as possi- 
ble, and using accessible volume techniques (Avni & Bah- 
call 1980) to produce unbiased estimators of the luminosity 
and dust mass functions. Currently we have completed one 
sample of 104 objects from the IRAS Bright Galaxy Sample 
(Soifer et al. 1989) selected at GO/xm. There is also an an 
optically selected sample (86 sources) taken from the CfA 



redshift survey (Huchra et al. 1983) for which data analysis 
is not yet complete (Dunne et al. 2001). This sample is de- 
signed to complement the IRAS sample by probing to lower 
submm luminosities, and also to address the possibility that 
the IRAS sample is missing a population of galaxies which 
are dominated by cold dust (and were therefore undetected 
by IRAS). 

The 850^m fluxes for the IRAS sample along with first 
estimates of the submm luminosity and dust mass functions 
have been presented in a previous paper (Dunne et al. 2000, 
henceforth Paper I). With the data then available to us 
(60, 100 and 850pim fluxes) we were only able to fit simple 
isothermal models of the dust spectral energy distributions 
(SEDs). While these provided a good empirical description 
of that data it was always possible that a more complicated 
multi-temperature model would be required in the presence 
of more FIR/submm data points. This would have an effect 
on the original dust masses we derived under the assumption 
of a single temperature, as at 850/im dust mass is roughly in- 
versely proportional to the assumed mean dust temperature. 
The idea of a multi-temperature model for dust in galaxies 
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is not new - it was initially invoked to explain relationships 
between IRAS observations and the HI, blue and prop- 
erties of the Milky Way and other galaxies (Cox, Krugel 
& Mezger 1986; de Jong & Brink 1987; Lonsdale Persson 
& Helou 1987; Boulanger & Perault 1988; Rowan-Robinson 
& Crawford 1989). In the most common version of these 
multi-temperature models, the FIR emission of a galaxy is 
hypothesised to originate in (i) a cool 'cirrus' component at 
T = 15 — 25 K, arising from diffuse dust associated with 
the Hi and heated by the general interstellar radiation field 
(ISRF) and (ii) a warm component at T > 30 K, from dust 
near to regions of high-mass star formation and Hll regions. 
The heating source for the warm component is Lyman alpha 
photons from young (OB) stars. As the longest passband on 
the IRAS satellite was 120/im and the cold 'cirrus' compo- 
nent is predicted to have a temperature of 15-25 K, only the 
warm component would have been detected by IRAS. Using 
IRAS 60 and lOO^m fluxes to determine the dust tempera- 
ture would therefore lead to an overestimate which, in turn, 
causes the dust masses to be underestimated (as at WOfxm 
Md oc T^^*~^^\ where /? is the dust emissivity index with a 
value believed to lie between 1 and 2). When the dust masses 
estimated using IRAS fluxes are combined with gas mass 
measurements, the gas-to-dust ratios are found to be a fac- 
tor of ~ 5-10 higher than for our own Galaxy (Devereux & 
Young 1990; Sanders et al. 1991), suggesting that IRAS does 
not account for all of the dust mass and strengthening the 
case for a two-component model. In order to demonstrate the 
presence of a cold component, measurements at wavelengths 
longer than 100/um are required (i.e. where the peak of the 
cold emission lies). However, even with a submm measure- 
ment at 800 — 1000/im, the dominance of the warm compo- 
nent in energy terms (L oc T'*^'') means that although most 
of the mass may reside in cooler dust, the excess emission 
produced by it at long wavelengths on the Ray leigh- Jeans 
tail is difficult to disentangle from a shallower wavelength 
dependence of the dust emissivity (i.e. a lower value of the 
dust emissivity index /3) Eales, Wynn- Williams & Duncan 
(1989); Chini et al. (1986). 

Due to the increase in spectral coverage brought by 
instruments such as COBE, ISO, SCUBA and the IRAM 
bolometer system, more recent studies have had access to 
measurements at multiple submm wavelengths in addition 
to IRAS fluxes. These have enabled astronomers to unam- 
biguously detect cold dust components at 15 < Td < 25 
K in several nearby spiral galaxies, including our own (e.g. 
NGC 891, NGC 5907, NGC 4656, M51, NGC 3627). Cold 
dust components have also been found in more IR lumi- 
nous/interacting systems such as IC 1623, NGC 4038/9, 
NGC 1068, NGC 3079 (Sodroski et al. 1997; Reach et al. 
1995; Ahon et al. 1998, 2001; Davics et al. 1999: Prayer et 
al. 1999; Papadopoulos & Seaquist 1999; Dumke et al. 1997; 
Braine et al. 1997; Haas et al. 2000; Neininger et al. 1996; 
Guolin et al. 1995; Sievcrs et al. 1994). The 850^m fluxes we 
presented in our first paper did not allow us to investigate a 
multi-temperature model because of the uncertainties in the 
dust emissivity index /3. We concluded that, in the absence 
of more submm data points, a single temperature SED with 
a /3 of 1.3 provided a good empirical description of the galax- 
ies in our sample. While we could not directly test whether 
there was a colder component present in our galaxies, we 
did discuss the implications - finding that our derived dust 



masses would increase by a factor ~ 2 and our gas-to-dust 
ratios would decrease to a value consistent with those deter- 
mined for the Galaxy and other spiral galaxies if there was 
a colder dust component present. The extra submm and 
mm data presented here allow a better determination of the 
SEDs of the galaxies in our sample than was possible in Pa- 
per I, and enable us to investigate the presence of colder 
dust in these bright IRAS galaxies. We use our improved 
knowledge of the dust properties of the galaxies to produce 
more reliable dust masses and a revised estimate of the dust 
mass function. Finally we compare the derived gas-to-dust 
ratios of the galaxies to the values for the Milky Way and 
other normal spirals. 

We wiU use Ho = 75 km s"^ Mpc"^ and go = 0.5 
throughout. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 Sample and Observations 

Paper I described our 850/im observations of a complete 
sample of 104 objects selected from the IRAS Bright Galgixy 
Sample (Soifer et al. 1989). We observed a subset of the BGS 
with SCUBA on the JCMT*, consisting of all galaxies with 
declination from —10° < 5 < 50° and with velocity > 1900 
km s~^, a limit imposed to try to ensure the galaxies fitted 
within the SCUBA field of view. The galaxies were observed 
between July 1997 and October 1998, using the jiggle-map 
mode of SCUBA. There are two arrays on SCUBA, 37 and 91 
bolometers for operation at long (850/im) and short (450/im) 
wavelengths respectively. They operate simultaneously with 
a field of view of ~ 2.3 arcmins (slightly smaller at 450/um). 
Typical beam sizes are ~ 15 arcsec and 8 arcsec at 850 
and 450/im respectively. Although the data at 450/im arc 
taken simultaneously with those at 850/im, they are only 
useful in dry and stable conditions {T450 < 1.5) meaning 
that only a small fraction (17 objects) of our sample has 
usable short wavelength data. We supplemented our 450/im 
data for the SLUGS galaxies with two observations from the 
JCMT axchive which we reduced, and also with submillime- 
tre, millimetre and FIR fluxes at A ~ 150 — 1400//m from 
the literature. 

2.2 Data Reduction and Analysis 

The method of data reduction and flux measurement has 
already been described in Paper I and also in more detail 
in Dunne (2000), but there were some slight differences at 
450/im which will be presented here. 

The data were reduced using the standard surf package 
(Jenness & Lightfoot 1998). The atmospheric extinction was 
measured using skydips which were performed at frequent 
intervals during observing (once every two or three hours). 
In order to estimate the zenith sky opacity, tx , a slab model 
for the atmosphere is fitted to the skydip measurements of 
the sky brightness temperature. A hot and cold load on the 

* The JCMT is operated by tlie Joint Astronomy Center on be- 
half of the UK Particle Physics and Astronomy Research Coun- 
cil, the Netherlands Organization for scientific Research and the 
Canadian National Research Council. 
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telescope are used to calibrate the skydip measurements, 
and the temperatures of these loads are used in the model 
fitting. The fits are generally good at 850/^m, but at 450/um 
a decent fit can be difficult to produce with the standard 
load temperatures. At the suggestion of Wayne Holland wc 
also performed fits with the following parameters altered: 
hot load temperature (Thot) changed from its default value 
to Tamb — 2 K, the ambient temperature, Tamb, being deter- 
mined from the FITS header of the data; cold load tempera- 
ture (Tcoid) changed to 95 K at 850/im and 170 K at 450/im. 
The combinations of default and new parameters provided 
four sets of model fits per skydip.^ Another estimation of r 
can be made from the Caltcch Submillimctrc Observatory 
(CSO) radiometer measurements at 225 GHz. The CSO ra- 
diometer performs a skydip every 10 minutes and relation- 
ships between this so-called Vcso' and the opacity values at 
450/850/im have been estabhshed by the JCMT staff. There 
is also a correlation between rgso and r45o which was used to 
obtain yet another estimate of the short wavelength opacity, 
as the 850/um skydips produced far more reliable fits than 
those at 450/im. All the estimates of tazo for each skydip 
were then averaged (after any obviously bad fits had been 
discarded) and that average value used in the extinction cor- 
rection. 

Fluxes were measured in apertures, chosen using the 

submillimetre and optical information as a guide and wher- 
ever possible matching that aperture used to determine the 
850/tm flux. The measurement uncertainties in the fluxes us- 
ing this method were described in Paper I, where we found 
through Monte-Carlo simulations that a straightforward ap- 
plication of 'CCD-type' shot noise greatly underestimated 
the true noise in the SCUBA maps. This is because the 
SCUBA maps are generally made with 1 arcscc pixels which 
do not represent independent regions on the sky. For a CCD 
the shot noise is given by tJahot = fpix-^/A^ap, whore fjpix is 
the standard deviation between pixels and iVap is the num- 
ber of pixels in the ob ject aperture. At 850/im we found that 
fshotlSSO) w 8 CTpix-^iVap. We performed the same Monte- 
Carlo procedure at 450/im and found that now 

<Tshot(450) W 4.4o-pix-\/ iVap 

The distinction between the shot noise corrections at 450 
and 850/tm is due to the differences in resolution at the two 
wavelengths; the ratio 4.4/8 is equivalent to the ratio of their 
beam sizes (8/15). The factors of 4.4 and 8 can be shown 
to be related to the size of the correlated regions within 
the maps as V A^corr , (A^'corr being the number of correlated 
pixels) so that Arcorr(850) w 64 and A"corr(450) w 19. These 
correlated regions are roughly 1/4 of the beam area, the 



^ Since all this data was reduced, the JCMT announced that 
prior to May 2000 the hot and cold load tcrnijcraturcs used by 
the skydip software had been inaccurate (Aiciiibald ct al. 2000) 
and that the standard value of T^ot had been too high. The rec- 
ommended changes should not greatly affect 850/im skydips but 
can make substantial differences to those at 450/tm depending on 
the conditions. As the suggested changes axe similar to what we 
did by varying the hot and cold temperatures, the r derived using 
the above method are not significantly affected by this realisation. 



reason for this lying in the linear weighting routines used by 
the rebinning task in SURF^ . 

Once the object flux had been measured in the aperture, 
it was calibrated in units of Jy using a map of a planet 
(Uranus or Mars) or the secondary submm calibrator, CRL 
618. The same aperture was used for the calibrator as for the 
object (paying attention to the orientation of the aperture 
relative to the chop throw, as the beam is elongated along 
the direction of the chop). The calibrated flux for the object 
is then given by 



Sobj(Jy) 



Sobj(V) 

Scal(V) 



X Scal(Jy) 



(1) 



where Scai(Jy) is the total flux of the calibration object in Jy 
which, for planets, was taken from the JCMT fluxes pro- 
gram. It was assumed that CRL 618 had fluxes of Sgso = 
4.56 Jy and S450 = 11-2 Jy (as listed on the JCMT cali- 
bration web page). Note that this is an inherently different 
method of calibration than that described in the mapping 
cookbook (Sandell 1997), which produces maps in units of 
Jy beam~^ by multiplying by a 'gain' or 'Flux Conversion 
Factor (FCF)'. This quantity (denoted C\) is simply the ra- 
tio of the flux of the planet in Jy beam~^ divided by the 
peaJi flux of the planet in Volts, as measured on the map. 



2.3 Calibration uncertainty at 450/im 

Even in the best atmospheric conditions, the zenith sky 
opacity of the atmosphere at 450/tm is ~ 0.6, and therefore 
the accuracy of the final 450/im fluxes depends critically on 
the accuracy of the extinction correction. In this section we 
describe a detailed investigation into the accuracy of 450/(m 
photometry with SCUBA. 



2.3.1 Extinction Correction 

The flux of an object at the top of the atmosphere is related 
to that received by the telescope by the following relation- 
ship 



-Atx 



where Foba,\ is the observed flux at the telescope, Fa,\ is the 
flux at the top of the atmosphere before attenuation, A is the 
airmass of the source and t\ is the zenith sky opacity. The 
opacity is wavelength dependent, being greater at shorter 
wavelengths and is a strong function of atmospheric water 
vapour content. Since r is directly related to the transmis- 
sion and thus to the relative flux observed, it is important to 
know its value as accurately as possible and to be measuring 
it often enough to keep track of any variability. At constant 
airmass, a change in t will produce the following change in 
the flux: 



El 
F2 



-Ari 



-At2 



(2) 



t For this reason the above relationships for shot noise (cshot) 
only apply to maps made in jiggle mode with 1 arcsec pixels and 
using the linear weighting routines 
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where At is defined as T2 — ti, and therefore Fi > F2. 
Expressed as a fractional change in flux, AF/F = (Fi — 
F2)/F\, this can be written as: 



25 



(3) 



We can make use of Eqn.3 to estimate tlic overall uncertainty 
in flux due to changes in t in the following way: An average 
value of r is calculated from the four skydip flts (varying 
the load temperatures) and the relationships with rg5o and 
Tcso if available. The standard deviation of these 5 or 6 esti- 
mates of T can be calculated, as well as the standard error 
on the mean value (cTrav)- This a^,^ can then be converted 
to a fractional error in the object flux using Eqn.3, where 
At — (7t,„. Fig. 1 shows AF/F plotted against r, using the 
450/im skydips at high r values and 850/im skydips for the 
lower T values. The uncertainty on the flux is a strong func- 
tion of T. On nights with high opacity the uncertainty in the 
fluxes is ~ 7—15%, and for t > 1.7 the inherent uncertainty 
on the flux due to the extinction correction makes observing 
a pointless exercise if accurate fluxes arc required. This can- 
not be overcome by making more frequent skydips (which is 
a separate issue), each point in Fig. 1 is for a single point in 
time. Fig. 1 therefore gives a lower limit to the extinction 
correction uncertainty. The problem of the sampling rate 
of T, i.e. how r varies between skydips, will be discussed 
later. In the era before skydips, the accepted practice was 
to observe the planetary calibrator at the same airmass as 
the target object. This helps to reduce the effects of uncer- 
tainties in T and would still be of benefit to observers using 
skydips today. One other point to note on observing tech- 
niques is that it is not necessarily wise to focus the telescope 
only immediately before planet measurements, as this means 
that the only maps in focus are those of the calibrators, not 
the targets. This will enhance any calibration problems from 
opacity correction and changes in the gain (see next section) . 



2.3.2 Changes m the gain 

The method used here to calibrate the fluxes relies on mea- 
suring the calibration object in the same aperture as the 
source. The variation of the calibration factor over time 
therefore sets an upper limit to the accuracy of the 450/Ltm 
galaxy fluxes. To asses this, all of the 450/Ltm calibration 
maps (13 maps from 11 nights) and a large number of 
850/im calibration maps (24 maps from 20 nights) were 
taken and the signal (in Volts) measured in different sized 
apertures (accounting for all the aperture sizes used to 
measure object fluxes). The 'aperture calibration factor', 
ACF — Seal (V) /Seal (Jy), was then calculated for each aper- 
ture on every map. The average ACF for each sized aperture 
at 450 and 850/im was taken and the standard deviation of 
all of the measurements used to give the uncertainty on that 
value. Both the average ACF and the uncertainties arc listed 
in Table 1 and it can be seen that the uncertainty is largest 
for the smallest apertures, but plateaus out at ~ 9 — 10 
percent at 450/im and 6 percent at 850/im. The last row 
(denoted Cx) lists the mean 'standard' flux conversion fac- 
tor (known as the FCF or gain) for the calibration maps. 
This is what is generally used to calibrate SCUBA maps, 
when the desired units are Jy beam~^. Table 1 shows that 
the values of Cx are far more variable than the aperture 
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Figure 1. Relative uncertainty on flux due to the uncertainty in 
deriving an accurate value for the sky opacity t, at a nominal 
airmass of 1.3. For observations at other airmasses (A) the y axis 
will scale as ~ ^4/1.3. The ordinate is a strong function of t and 
places a limit on the sky conditions which are usable for 450/im 
observations when accurate calibration is required. The crosses 
represent 450/im r values, the stars are 850/tm values. The scatter 
and relative uncertainty at 850/im are much lower than at 450/im. 
The fitted Une has the form y = -0.58 + 5.24t - 6.97t'^ + 5.56t^. 



calibration factors (ACF) used here to convert Volts to Jy. 
This is because measuring Cx uses the peak flux, which will 
be strongly affected by changes in beam shape due to sky 
noise, pointing drifts, chop throw and the dish shape (most 
critical at 450/im). Anything which broadens the beam will 
remove flux from the peak and place it further away, mak- 
ing Cx higher. The inherent surface accuracy and thermal 
fluctuations in the shape of the dish cause large variations 
in C'x at 450/im (the changes arc still present, but smaller 
at 850/tm); thus the gain at 450/tm often does not stabilise 
until after midnight, and can vary by more than a factor of 
two in the course of one night. 

The results in Table 1 assume that long term changes 
in the instrumental response and manual alterations to the 
dish shape (via holography), are smaller than the short 
term changes in the beam pattern due to thermal fluctu- 
ations and other effects. The validity of this assumption was 
investigated by comparing observations made within ~ 1 
month of each other to see if the variation was any less. At 
450/tm this was not the case and so the values from Ta- 
ble 1 represent the basic uncertainty in the ACF. However, 
at 850/im, a significant reduction was found when using 
observations taken closer together (Table 2). The changes 
in Cx over the 14 month period tie in with known im- 
provements/deteriorations in the dish surface. Therefore, at 
850/tm the values from Table 2 should be used to estimate 
the uncertainty in the ACF. 

The ACF method presented here provides a more ro- 
bust way to calibrate fluxes, which is consistent to within 
10 percent at 450/im and 2.6 percent at 850/im (in a typical 
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Table 1. Relative uncertainties in gains and calibration factors used in flux measurement 



5 



Aperture SSO/im (24) 450/im (10) Notes 





mean 


a 


% a 


mean 


a 


% a 




15" 


2.229 


0.167 


7.5 


11.13 


1.37 


12 


Aperture 


30" 


1.134 


0.080 


7.0 


7.20 


0.63 


9 


calibration 


45" 


1.013 


0.060 


6.0 


5.79 


0.49 


9 


method 


60" 


0.942 


0.052 


5.5 


5.12 


0.50 


10 




90" 


0.879 


0.050 


5.6 


4.75 


0.50 


10 




Ca 


278 


22 


8.0 


815 


144 


18 


'Standard' peak gain or FCF 



Mean, standard deviation and % variation in the 'aperture calibration factors' and the 'standard' peak calibration factor for both 
wavelengths over the entire course of the observations. Numbers in parentheses are the numbers of maps used to determine these 
figures. The method of calibration is in the final column. At 450/im, the variation here is no greater than that over a shorter time 
period but for 850//m it is an overestimate, see Table 2. 

Table 2. Variation in gains and calibration factors at 850/im over shorter time scales than Table 1. 



Dates Nmaps C\ % change in ACF 







Mean 


% a 


15" 


30" 


45" 


60" 


90" 


July/ Aug 97 


4 


266 


3.0 


2.5 


2.3 


1.4 


0.4 


0.5 


Dec/ Jan 97 


4 


283 


4.0 


5.0 


4.4 


4.4 


3.6 


4.1 


March 98 


9 


259 


4.3 


3.1 


2.1 


1.8 


1.8 


2.3 


June/July 98 


5 


293 


6.0 


4.5 


3.8 


2.9 


2.3 


2.4 


Sept 98 


2 


294 


5.7 


5.2 


4.1 


2.8 


1.7 


1.3 


Average 




4.6 


4.0 


3.3 


2.6 


2.0 


2.1 



850/im % variation in C\ and ACF for time periods of ~ 1 month. The mean C\ can be seen to change on these time scales due to 
deliberate alterations in the dish shape. As a result of this, the variations within these time frames are smaller than those seen in 
Table 1 and thus is it the values from this Table which should be used in estimating (Jgain at 850/im. 

eludes the arehival data). Fluxes from the literature for the 
SLUGS gala^xies are listed in Table 4 and are from a variety 
telescopes (ISO, JCMT, IRAM, OVRO). Other fluxes for 
Arp220 are listed separately in Table 5. There are complica- 
tions when making comparisons between fluxes measured at 
different wavelengths if there are significant beam size dif- 
ferences between them. We did not include fluxes from the 
literature if the beam sizes were smaller than the extent of 
the submm emission. 

3.1 Corrections applied to fluxes 

3.1.1 SCUBA maps 

We measured the 450/im emission in the same sized aper- 
ture used for the 850//m observations, which was chosen 
to maximise the signal-to-noise (S/N) at the longer wave- 
length. However, because of the lower S/N at the shorter 
wavelengths, the 450/im emission is sometimes only clearly 
detected in the central region of the aperture, and system- 
atic effects, such as those caused by poor sky subtraction 
are more of a concern. We investigated this by convolving 
the 450 and 850/tm maps to the same resolution and tak- 
ing a small aperture around the bright central region on 
both maps. The ratio of the 450/tm and 850/tm fluxes in 
this region was compared to that measured using the larger 
aperture. If the ratio for the larger aperture was significantly 



aperture), over a time period of ~ 1 month at 850/iin and 
much longer at 450/im, and in a wide range of conditions. 

The total calibration uncertainty at 450/tm can be es- 
timated as follows: 

/ 2 , 2 ,2 ,2 n1/2 / .\ 

fcal = (Cr + fgain + 0"vary + f"abs j (,4) 

where <Jt- is the uncertainty in r from skydip fitting (on av- 
erage 7 per cent for these observations, sec Fig. 1), o-gain is 
the uncertainty in the ACF factor (on average 9 per cent, sec 
Table 1), (Jabs is the absolute uncertainty in the brightness 
temperatures of the planets which are used for calibration 
(estimated to be 5 per cent (SUN/213)) and avary repre- 
sents the uncertainty in how the sky opacity has changed 
between skydips. This is not really calculable, as it will be 
a function of skydip frequency and atmospheric instability. 
However, the calibrators themselves have had their fluxes 
corrected for extinction, and so (Tvary must be part of again. 
We estimate an upper limit on o"vary as a^/ ""gain ~ '^^ which 
is ~ 2 — 4 percent. Including all terms, the total calibra- 
tion uncertainty is ~ 10 percent at 850/tm and 15 percent 
at 450/im. 



3 RESULTS 

Table 3 lists the 450/im fluxes which we analysed and which 
were measured in the way described in Section 2.2 (this in- 
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Table 5. Other FIR/submm fluxes for Arp220. 



Wavelength Flux Ref. 
(/im) (Jy) 



60 109.4 g 

65 114.4 g 

90 95.2 g 

100 105.3 g 

120 73.4 g 




150 


58.5 


g 


170 


47.9 


g 


180 


37.5 


g 


200 


34.8 


g 


350 


11.7 


e 


450 


6.286 


* 


850 


0.832 


* 


1100 


0.350 


e 


1250 


0.226 


c 


1300 


0.175 


h 


1400 


0.140 


1 



References as for Table 4, * submm fluxes from this work. ISO 
fluxes (g) are colour corrected as described in that paper. 

lower than the central ratio, the total 450/um flux would be 
corrected by a factor which would reproduce the central ra- 
tio across the whole galaxy. In doing this, it is implicitly as- 
sumed that the true flux ratio is constant across the gala^xy. 
Observations of the large, well-resolved galaxies NGC 891 
and NGC 7331 (Alton et al. 1998; Bianchi et al. 1998; Israel 
et al. 1999; Alton et al. 2001) and of the interacting pair 
VV114 (Fraycr ct al. 1999), and of the Milky Way (Sodroski 
et al. 1997) suggest that the actual change in S450 /S'sso due 
to dust temperature gradients would be only ~ 14 percent 
or less. In practise, such a correction was only necessary for 
four objects, two of which (NGC 5962 and NGC 7541) were 
large and extended over most of the array. The other two 
(UGC 2369 and UGC 2403) were observed in quite noisy sky 
conditions on the same night. The corrections were not very 
large, ranging from 10-34 per cent and they do not affect 
the conclusions. Fluxes in Table 3 have been corrected for 
this eflfect, and any corrections used are also listed. 



3.1.2 Literature fluxes 

Due to the uncertainty involved in making beam cor- 
rections, only literature fluxes which were obtained with 
single-element detectors with beams larger than the area of 
SCUBA emission (i.e. long wavelength ISO data, observa- 
tions of very small sources) or which come from maps were 
used. However, a small correction was applied to the fluxes 
for NGC 4418, which was observed at 450, 800 and llOO^tm 
by Roche & Chandler (1993) using UKT14 on the JCMT. 
This source is very compact with almost all of the flux in 
one SCUBA beam. The SCUBA map was used to find the 
850/xm flux in the equivalent region to the UKT14 beam at 
each wavelength. In this way correction factors of 1.10, 1.25 
and 1.09 were derived for the 450, 800 and llOO/nm fluxes 
respectively. This is a small difference and has very little 
effect on the derived parameters in later discussions. 



(a) NGC 958: cold (b) UGC 2982: warm 

Figure 2. a) SED of galaxy with a prominent cold component, 
not compatible with a single temperature SED. b) SED of galaxy 
with dominant warm component, any cold dust is very difficult 
to see in the SED (if it is there at all). The solid lines indicate the 
composite two-component SED (parameters listed in Table 6). 
The dot-dashed lines represent the warm and cold components. 



3.2 Two-component Spectral Energy 
Distributions. 

The emission at a particular frequency will now be repre- 
sented as the sum of two modified Planck functions, each 
with a different characteristic temperature. In reality, the 
dust is probably at a range of temperatures reflecting the lo- 
calised heating sources in the ISM, but the approximation of 
characteristic temperatures for the warmer and cooler dust 
is not a bad one. For the optically thin regime, this can be 
expressed as: 

S, = N^x i^'^B{i^, T^) -I- iVc X z/''B(zy, Te) (5) 

where A'^^, and Nc represent the relative masses in the warm 
and cold components, T-u, and Tc arc the temperatures, 
B(i/, T) is the Planck function for each component, and /3 
is the dust emissivity index (assumed to be the same for 
each). This model was fitted to the fiuxes and the parame- 
ters which produced the minimum were found. Initially 
the value of /3 was allowed to vary between 1 — > 2; r,„ was 
constrained by the IRAS 25pm flux, which was not allowed 
to be exceeded in a fit (N.B. a fit was not made to this point 
as this would require an additional hot component); Tc was 
allowed to take any value lower than T^,. For sources with 
only 4 fluxes (not including 25/jm), there are not enough 
data points to provide a well-constrained fit, consequently 
the values of Xmin s^re unrealistically low and a large range 
of parameters will provide very acceptable fits. The param- 
eters producing the best fits are in Table 6 along with the 
original isothermal SED parameters from Paper I - derived 
from fits to the 60, 100 and 850/xm points only. Uncertain- 
ties in the parameters from the two-component fits are also 
given and these were produced using the bootstrap method, 
as described in Paper I. 

The relative contribution of the cold component to the 
SEDs shows a large variation (the best parameter to describe 
this would be Nc/N,^). Fig. 2a shows an example of a very 
'cold' SED where the cold component is clearly visible and, 
once the 450/Lim point is included, a single temperature fit 
is excluded. This is in contrast to other objects (Fig. 2b) 
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Table 3. Flux densities for the SLUGS galaxies 



(1) 


(2) 


(3) 




(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Name 


R.A. 


Dccl 




C2 




5100 


S'45() 


(^450 


Ssso 


O-850 


fcor 


^850 




(J2000) 


(J2000) 


(km s~^) 


(Jy) 


(Jy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 






UGC 903 


1 21 47.9 


+17 


35 


34 


2518 


7.91 


14.58 


1500 


315 


178 


26 




8.43 


NOG 958 


2 30 42.8 


—02 


56 


23 


5738 


5.90 


14.99 


2251 


428 


262 


34 




8.60 


UGG 2369 


2 54 01.8 


+14 


58 


14 


9400 


7.68 


11.10 


523 


120 


72 


13 


1.11 


7.26 


UGC 2403 


2 55 57.2 


+00 


41 


33 


4161 


7.51 


11.77 


1010 


202 


111 


18 


1.26 


9.10 


NGC 1614^ 


4 34 00.0 


—08 


34 


45 


4778 


33.12 


36.19 


981 


167 


140 


20 




7.45 


NGC 1667 


4 48 37.2 


-06 


19 


12 


4547 


6.24 


16.54 


1183 


272 


163 


22 




7.26 


NGC 2856 


9 24 16.2 


+49 


14 


58 


2638 


6.15 


10.28 


993 


268 


89 


16 




11.16 


NGC 2990 


9 46 17.2 


+05 


42 


33 


3088 


5.49 


10.16 


1275 


332 


110 


19 




11.60 


UGC 5376 


10 00 26.8 


+03 


22 


26 


2050 


5.94 


11.49 


1258 


415 


148 


23 




8.50 


ARP 148 


11 03 54.0 


+40 


50 


59 


10350 


6.95 


10.99 


646 


156 


92 


20 




7.02 


MCG+00-29-023 


11 21 12.2 


-02 


59 


03 


7464 


5.40 


8.87 


571 


166 


84 


13 




6.80 


ZW 247.020 


14 19 43.3 


+49 


14 


12 


7666 


5.91 


8.25 


284 


111 


36 


8 




7.90 


1 Zw 107^ 


15 18 06.1 


+42 


44 


45 


11946 


9.15 


10.04 


423 


93 


60 


14 




7.05 


IR 1525+36 


15 26 59.4 


+35 


58 


37 


16009 


7.20 


5.78 


252 


70 


33 


8 




7.64 


ARP 220 


15 34 57.2 


+23 


30 


11 


5452 


103.33 


113.95 


6286 


786 


832 


86 




7.56 


NGC 5962 


15 36 32.0 


+16 


36 


22 


1963 


8.99 


20.79 


1959 


372 


317 


37 


1.12 


6.18 


NGC 6052 


16 05 13.0 


+20 


32 


34 


4712 


6.46 


10.18 


721 


230 


95 


15 




7.59 


NGC 6181 


16 32 21.2 


+19 


49 


30 


2379 


9.35 


21.00 


1470 


456 


228 


37 




6.45 


NGC 7541 


23 14 43.4 


+04 


32 


04 


2665 


20.59 


40.63 


2639 


686 


427 


60 


1.34 


6.18 



(1) Most commonly used name taken from the IRAS BGS (Soifer et al. 1989). (2) Right ascension J2000 epoch. (3) Declination J2000 
epoch. (4) Recession velocity taken from NED. [The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion 

Laboratory, California Institute of Tcclmology, under contract with the National Aeronautics and Space Administration.] (5) 60/im flux 
from BGS. (6) lOO^m flux from BGS. (7) 450^m flux (this work). (8) error on 450/im flux, calculated in the manner described in 
Section 2.2 and inclusive of a 15% calibration uncertainty. (9) 850/im flux (Paper I). N.B. The values for NGC 1614 and IR 1525+36 
have been revised from those listed in Paper I in light of new data. (10) error on 850/im flux (Paper I, includes a 10% calibration 
uncertainty). (11) Correction factor applied to 5450 if a decline in S'45o/S'850 was found (Section 3.1); corrected values arc used in Cols 
5,7,12 and later in this paper. (12) Ratio of 450(Um to 850//m fluxes. ^ indicates that the data was taken from the JCMT archive.. 



Table 4. Other submillimetre fluxes for the IRAS sample from the literature. 



(1) 


(2) 


(3) 




(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Name 


R.A. 


Dccl 




cz 




Sim 


5850 




A 


Sx 


Refs. 




(J2000) 


(J2000) 


(km s-l) 


(Jy) 


(Jy) 


(mJy) 


(mJy) 


{urn) 


(Jy) 




NGC 520 


1 24 34.9 


+03 47 


31 


2281 


31.55 


46.56 


325 


50 


1200 


0.076 


a 


UGC 2982 


4 12 22.5 


+05 32 


51 


5305 


8.70 


17.32 


176 


34 


180, 190, 1250 


9.9, 9.2, 0.047 


b,c 


NGC 2623 


8 38 24.1 


+25 45 


16 


5535 


25.72 


27.36 


91 


14 


350, 750 


2.225, 0.170 


d 


NGC 3110 


10 04 02.0 


-06 28 


31 


5048 


11.68 


23.16 


188 


28 


1250 


0.048 


c 


IR 1017+08 


10 20 00.2 


+08 13 


34 


14390 


6.08 


5.97 


36 


6 


350, 750 


0.450, 0.049 


d 


IR 1056+24 


10 59 18.2 


+24 32 


34 


12912 


12.53 


16.06 


61 


13 


350, 450, 750 


1.24, 0.533, 0.085 


d 


IR 1211+03 


12 13 46.1 


+02 48 


40 


21703 


8.39 


9.10 


49 


10 


450 


0.429 


d,e 


NGC 4418* 


12 26 54.7 


-00 52 


39 


2179 


42.32 


30.76 


255 


37 


450, 800, 1100 


1.474, 0.300, 0.093 


f 


UGC 8387 


13 20 35.3 


+34 08 


22 


7000 


13.69 


24.90 


113 


15 


350, 750 


1.946, 0.196 


d 


Zw 049.057 


15 13 13.1 


+07 13 


31 


3927 


21.06 


29.88 


200 


27 


1250 


0.047 


c 


NGC 7592 


23 18 22.1 


-04 24 


58 


7350 


8.02 


10.50 


108 


19 


80, 180, 200 


11.2, 5.8, 4.6 


b 


NGC 7679 


23 28 46.7 


+03 30 


41 


5138 


7.28 


10.65 


93 


15 


80, 180, 200 


9.4, 5.8, 4.6 


b 


NGC 7714 


23 36 14.1 


+02 09 


18 


2798 


10.52 


11.66 


72 


13 


80, 180, 200 


11.0, 4.9, 4.3 


b 



(1) Most commonly used name taken from the IRAS BGS (Soifer et al. 1989). (2) Right ascension J2000 epoch. (3) Declination J2000 
epoch. (4) Recession velocity taken from NED. (5) 60(Um flux from BGS. (6) WOfMa flux from BGS. (7) 850/xm flux (Paper I). (8) error 
on 850^m flux (Paper I, includes a 10% calibration uncertainty). (9) Wavelengths of measurements (microns). (10) Fluxes at each A in 
Jy. (11) References for fluxes: a) Braine & Dumke 1998; b) Sicbenmorgen et al. 1999; c) Carico et al. 1992; d) Fox 2000; e) Rigopoulou 
et al. 1996; f) Roche & Chandler 1993; g) Klaas et al. 1997; h) Downes & Solomon 1998; i) Woody et al. 1989 

* Fluxes for NGC 4418 have been beam corrected using the SCUBA 850/xm map (Section 3.1). Corrections were small (10% increase at 
450jum, 25% at 800//m and 9% at 1100). 
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Table 6. SED parameters for isothermal and free /3 two-component fits. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Name 




Td 


02 








Tc 






tn 


Nx 






(K) 






(K) 


(K) 


(K) 


(K) 








UGC 903 


1.2 


34.4 


2.0 


(0.36) 


45 


(5.1) 


21 


(4.2) 


91 


(52) 


4 


NGC 958 


1.2 


30.8 


2.0 


(0.25) 


44 


(9.1) 


20 


(2.9) 


186 


(341) 


4 


UGC 2369 


1.4 


36.2 


1.5 


(0.35) 


41 


(13.3) 


32 


(6.0) 


3.3 


(101) 


4 


UGC 2403 


1.2 


36.8 


2.0 


(0.30) 


50 


(6.5) 


22 


(4.8) 


100 


(65) 


4 


NGC 1614 


1.6 


39.2 


1.6 


(0.23) 


40 


(9.6) 


31 


(9.1) 


0.06 


(20) 


4 


NGC 1667 


1.6 


28.4 


1.6 


(0.27) 


31 


(10.2) 


28 


(4.4) 


3.9 


(136) 


4 


NGC 2856 


1.3 


35.0 


2.0 


(0.27) 


45 


(8.7) 


22 


(4.5) 


58 


(134) 


4 


NGC 2990 


1.3 


33.8 


2.0 


(0.17) 


42 


(8.9) 


21 


(2.6) 


58 


(133) 


4 


UGC 5376 


1.2 


33.8 


2.0 


(0.42) 


44 


(7.4) 


21 


(5.5) 


87 


(116) 


4 


Arp 148 


1.3 


35.6 


1.5 


(0.35) 


47 


(7 P.) 


30 


ffi 5~) 

(V.O) 


15.5 


CSTl 
Voij 


4 


MCG+00-29-023 


1.3 


35.0 


1.5 


(0.32) 


40 


(9.3) 


27 


(5.8) 


4.8 


(50) 


4 


Zw 247.020 


1.7 


34.4 


1.7 


(0.24) 


36 


(14.0) 


28 


(5.4) 


0.4 


(97) 


4 


1 Zw 107 


1.3 


41.6 


1.4 


(0.39) 


43 


(10.5) 


27 


(8.3) 


0.9 


(30) 


4 


IR 1525+36 


1.1 


53.0 


1.7 


(0.39) 


57 


(5.8) 


26 


(9.7) 


15.2 


(17) 


4 


NGC 5962 


1.2 


32.0 


1.4 


(0.28) 


32 


(9.8) 


21 


(7.3) 


1.0 


(77) 


4 


NGC 6052 


1.3 


35.6 


1.7 


(0.34) 


55 


(12.8) 


26 


(6.8) 


65 


(65) 


4 


NGC 6181 


1.4 


30.8 


1.4 


(0.33) 


32 


(11.6) 


31 


(5.3) 


1.0 


(165) 


4 


NGC 7541 


1.3 


33.2 


1.4 


(0.31) 


34 


(9.6) 


27 


(4.8) 


0.6 


(63) 


4 


NGC 520 


1.4 


36.2 


2.0 


(0.24) 


50 


(6.4) 


24 


(4.4) 


61 


(42) 


4 


UGC 2982 


1.4 


32.0 


1.6 


(0.27) 


33 


(10.7) 


22 


(6.7) 


1.6 


(872) 


6 


NGC 2623 


1.6 


39.8 


2.0 


(0.15) 


50 


(5.7) 


27 


(1.5) 


18 


(17) 


5 


NGC 3110 


1.5 


32.0 


1.9 


(0.30) 


31 


(12.4) 


18 


(10.1) 


2.7 


(83) 


4 


IR 1017+08 


1.2 


45.2 


1.3 


(0.32) 


44 


(8.2) 


27 


(13.3) 


0.3 


(18) 


5 


IK, lUoD+z4 


1.7 


35.6 


z.O 


(0.23) 


40 


/I o 1 ^ 
(12.1) 


26 


(b.b) 


5.6 


(67) 


6 


IR 1211+03 


1.3 


42.2 


1.9 


(0.34) 


49 


(5.5) 


25 


(7.5) 


20 


(44) 


4 


NGC 4418 


0.8 


62.0 


1.5 


(0.35) 


55 


(6.3) 


23 


(6.2) 


8.9 


(25) 


6 


UGC 8387 


1.8 


30.8 


1.8 


(0.16) 


32 


(11.3) 


31 


(3.8) 


0.7 


(1340) 


5 


Zw 049.057 


1.4 


36.2 


2.0 


(0.33) 


44 


(5.0) 


23 


(7.0) 


28 


(35) 


4 


Arp 220 


1.2 


42.2 


1.8 


(0.24) 


47 


(3.6) 


19 


(3.6) 


19 


(14) 


16 


NGC 7592 


1.1 


39.8 


1.5 


(0.33) 


39 


(11.0) 


20 


(11.1) 


3.5 


(76) 


6 


NGC 7679 


1.2 


38.0 


1.5 


(0.29) 


48 


(11.8) 


29 


(6.5) 


13.7 


(64) 


6 


NGC 7714 


1.3 


41.0 


1.5 


(0.30) 


60 


(12.8) 


32 


(7.9) 


19.0 


(171) 


6 



(1) Galaxy name. (2) Emissivity index derived from an isothermal fit to tlic 60, 100 and 850/im fluxes (Paper I). (3) Dust temperature 
derived from an isothermal fit to the 60, 100 and 850/im fluxes (Paper I). (4-5) Best-fitting value of dust emissivity index j3 using a 
two-component SED and la uncertainty from bootstrap calculation. (6-9) Fitted warm and cold component temperatures with lu 
uncertainties. (10-11) The ratio of cold dust mass to warm dust mass and la uncertainty. (12) Number of fluxes used in the 
two-component fit (not including the 25/Lfm point). 

Galaxies below the line are those with extra FIR/submm fluxes taJcen from the literature while above are those for which we have our 
own data. 

where a cold component can be included but its contribution 

is barely noticeable and it is statistically indistinguishable 
from an isothermal SED. The problem with indeterminate 
objects such as that in Fig. 2b is not simply due a limited 
number of data points. For example, the SED of Arp 220 
(Fig. 3) has more than enough data points, yet can still 
be fitted by either single temperature models with lower 
(3 or two-component models with higher /3 (however two- 
component models produce the best fits). 

The distribution of /3 for the sample is presented in 
Fig. 4, where the original single temperature f)s for these 
galaxies are also shown. The two-component P distribution 
(solid) is shifted toward higher values than the single tem- 
perature one (dot-dashed) and also seems to have two mgix- 
ima, one at /3 = 2 and another at /? = 1.5. The Tc distribu- 
tion (Fig. 5) is also quite odd with 2-3 peaks at ~ 22, 27 and 



31 K. The galaxies producing the /? w 1.5 peak seem to be 

the same ones responsible for the higher Tc peaks at 27-31 
K. This bi-modal behaviour shows no correlation with other 
properties such as Lflr- 

3.3 Is there a Cold Component? 

To try to determine whether or not a cold component is nec- 
essary, we also tried to fit isothermal SEDs to the expanded 
data sets for various assumptions about (3. Table 7 shows 
how many good fits were produced for freely varying /3 as 
well as for fixed values of 1.5 and 2.0. It can be seen that 
if no constraints arc placed on /3 then only 7 objects are 
incompatible with a single temperature SED and require a 
cold component (although many more will have better fits 
to a two-component model than to a single temperature one. 
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Wovelength (urn; 

Figure 3. Best fitting SEDs for Arp 220 wiicn /3 is a free param- 
eter. Although there arc more than enough data points determine 
the free parameters Ty,, Tc, (3 and N^/Nw it is still possible to fit 
SEDs with either single temperatures (dashed) and lower /3, or 
two temperatures and steeper j3 (solid). The fundamental prob- 
lem in confirming the existence of colder components is not sim- 
ply one of a lack of data points, but rather a degeneracy in the 
parameters. The parameters for this fit are listed in Table 6. 



1.5 

emissivity index {Ji) 



Figure 4. Distribution of /3 for the 32 sources with subiriiii fluxes 
at more tliaii ouc wavelength. The solid line is the 2 component 
/3 distribution and the dot-dash line represents the /3 from the 
original single temperature fits to the 60, 100 and 850/im fluxes. 



cold temperature (K) 

Figure 5. Distribution of cold temperatures when /? is allowed 
to take on values between 1-2. 



Table 7. Galaxies which could be fitted by single temperature 
SEDs for various /3. 



p 


Fit (xl < 1) 


No fit 


free 


25 


7 


1.5 


18 


14 


2.0 


2 


30 



3.4 Constraints on (5 

The value of fH is of great importance, as it tells us about 
the behaviour of the dust emissivity with wavelength. The 
dust mass opacity coefficient Kd has only been mesisured at 
FIR wavelengths (~ 125/im) and must be extrapolated to 
the wavelength of observation as oc \~^. Knowledge of /? 
therefore improves our estimate of Kd(850), and hence dust 
mass. A well determined l3 is also necessary for interpreta- 
tion of the high redshift submm observations, where only 1 
or 2 fluxes are available, making a simultaneous determina- 
tion of T and (3 impossible and attempts to use the submm 
flux ratios as redshift indicators very uncertain. However, 
one must also be wary that there may be evolution in the 
dust properties, such as temperature, emissivity and extinc- 
tion behaviour with redshift and metallicity. These issues are 
not yet fully understood, but it is worth noting that obser- 
vations of high-redshift radio galaxies and quasars are ade- 
quately represented by an emissivity proportional to j/i-5-2 
(Benford et al. 1999; Ivison et al. 1998; Downes et al. 1992; 
Hughes et al. 1997; Chini & Kriigel 1994). Section 3.2 has 
shown that two-component fits to small data sets with a free 
/3 cannot tell us anything conclusive about either the value 
of (3 or the putative cold compoiicut. We will now explore 
other ways to deduce the value of (3. 



it is not possible to rule the single temperature model out). 
The balance shifts as 3 is forced to be steeper, with only 
two objects allowing single temperatures for (3 = 2. This 
suggests that if we knew the value of (3 (if indeed one value 
should be applied to all of the sources), we could determine 
whether or not a cold dust component is necessary for these 
galaxies. 



Theory and observations 

A mixture of laboratory experiments and theoretical models 
have produced various estimates of what should be for dif- 
ferent chemical compositions believed to represent interstel- 
lar dust, at various temperatures and sizes. These studies are 
not in great consensus. Draine & Lee (1984) measured opti- 
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cal constants for a mixture of silicates and graphite and in 
their model /3 ~ 2 is applicable for 40 < A < lOOO/xm (Gor- 
don 1995). Agladze et al. (1996) found in laboratory experi- 
ments that there may be a dependence of /3 on Td for certain 
types of amorphous silicates; /3 = 2 at 20 K and decreases 
with temperature although it is still > 1.5 at 30 K (however 
these measurements were made in the millimetre regime not 
the submillimetre). Other samples in the same experiment 
showed no such behaviour with temperature, and had a con- 
stant value of B — 1.2. Mcnnclla ct al. (1998) also found a 
dependence of /3 on T between 20/Lim - 2 mm but at higher 
temperatures (24 - 295 K). They used samples of silicates 
and amorphous carbons believed to be analogous to inter- 
stellar dust. The values of /3 appropriate for the silicates 
in the temperature range appropriate for ISM dust grains 
(10-50 K) are still clustered around 2 without a very strong 
temperature trend. For the amorphous carbon grains a low 
/? ~ 1 was favoured, while graphitic grains had /3 = 2 (Men- 
nella 1995), however the low /3 for the amorphous carbons 
could be due to the very small particle sizes used < O.Qljj,m. 
Fluffy composite or fractal grains have also been consid- 
ered, which could be representative of grains found in dark 
clouds. Mathis & Whiffen (1989) found (3 = 1.5 for compos- 
ite grains and Wright (1987) gives 0.6 < l3 < 1.4 for fractal 
grains. However, similar fluffy composites were investigated 
by Kriigel & Siebenmorgen (1994) who found /3 = 2, except 
for the largest grains (> 30/im) which are not believed to 
be representative of general diffuse ISM dust anyway. From 
these studies, it is difficult to justify a universal value of (3 
for diflerent regions in the ISM, since composition, size and 
possibly temperature could have an effect. 

Observations of diverse environments within the Galaxy 
also suggest a range of values for /3. COBE/FIRAS studies 
of the diffuse ISM in the Milky Way (Sodroski et al. 1997; 
Reach et al. 1995; Masi et al. 1995) find /? = 1.5 - 2.0 de- 
pending on whether one or two temperature components are 
fitted (a two-component model with /3 = 2 gives the best fit). 
Warm clouds associated with Hll regions are also found to 
have (3 = 2 (Gordon 1988). Recent results from the submm 
balloon experiment PRONAOS (Dupac et al. 2001) find ten- 
tative evidence for an inverse dependence of 13 on T but the 
result is largely based on one cloud (OMC-1) having a lower 
(3 and higher temperature than the others. There axe prob- 
lems with this interpretation if OMC-1 contains dust at more 
than one temperature (i.e. the same problem we have here), 
also the temperatures of the clouds investigated by Gordon 
(1988) ranged from 30-80 K, yet all apart from OMG-1 had 
/3 = 2. It may be that OMC-1 contains a grain population 
with larger than normal sizes, or more than usual amounts 
of a particular grain composition. Colder clouds in Orion A 
have also been observed, and /3 found to be 1.9 (Ristorcelli et 
al. 1998). However, observations of both young and evolved 
stars have revealed low (3 ~ 0.2— 1.4 (Weintraub et al. 1989; 
Knapp et al. 1993) and this is usually attributed to grains 
growing to very large sizes in cool, dusty envelopes or disks. 

In contrast, multi- wavelength observations of external 
galaxies are almost universal in their finding that 1.5 < /3 < 
2 and favouring (3 = 2 (Chini et al. 1989; Chini & Kriigel 
1993; Alton et al. 1998; Frayer et al. 1999; Bianchi et al. 
1998; Braine et al. 1997). It may be that the contribution to 
the global SED from dust around stars (which displays low 
(3) is small when averaged over the whole galaxy. Assuming 



that the production of dust occurs in similar ways in all 
galaxies (analogous to assuming star formation mechanisms 
are universal) then maybe we should expect the global dust 
properties to be generally consistent from galaxy to galaxy 
(after all we must make some such assumption if we wish to 
compare dust masses etc). 

Colour plots: the ratio 84^,0 /Ss5o 

Figure 6 shows the colour-colour plot for 5'6o/<S'450 versus 
Seo/Ssso for the objects in this sample which have 450/im 
fluxes (stars and circles) plus some objects from the litera- 
ture (which are not in our sample) with fluxes at these wave- 
lengths (IC 1623, NGC 891, Mrk 231, Mrk 273, UGC 5101 
- open triangles) plus some of the optically selected galaxies 
from our as yet incomplete sample (diamonds). The stars 
refer to objects which we observed and reduced, and the 
circles are 450/im fluxes from the literature for this sample. 
There arc a few points to note here: 

(i) The relationship is very tight (correlation coefficient 
rs = 0.96, significance=5.8a) and the scatter is consistent 
with the uncertainties in the flux ratios (denoted on the 
figure). 

(ii) This relationship holds across 2-3 orders of magni- 
tude in Lfir, even though there are many different types of 
galaxies, ranging from optically selected galaxies with low 
Lfir to very active ULIRGS and AGN such as Arp 220 and 
Mrk 231. 

(iii) The best fitting line (least squares) is: 

log(S'6o/S'45o) = (1.01±0.03) log(S'6o/S'8Bo)-(0.909±0.045)(6) 

This can be re-written as S'6o/<S'450 = 0.123 (Seo/S'sso)^"^ 
which implies that, within the uncertainties, the ra- 
tio 5450 /5850 is constant. From the fluxes, the mean 
5'45o/S'86o = 7.9 with a = 1.6 which is a remarkably tight 
distribution, consistent with the errors on the fluxes. 

What does a constant value of S450/S850 mean? The 
flux ratio can be expressed as a function of T and /3 as 
follows: 

^9450 _ / t^450 \ ^ ^ 6(^-450, r) 

Since the submm fluxes are on the Ray leigh- Jeans side of 
the Planck function the temperature dependence is not very 
steep, roughly oc T. A constant value for the ratio therefore 
implies a small range of real (3 and Tc in the ISM of galaxies. 
We further investigated the implications of the relationship 
in Fig. 6 in the following way: Hypothetical SEDs were cre- 
ated using the parameters T™, Tc, Nc/N^ and /3, the values 
for these being either specified, drawn from a uniform ran- 
dom distribution between two given limits, or drawn from 
a Gaussian distribution with a designated mean and disper- 
sion. The various models are listed in Table 8. Each SED 
was used to produce fluxes at 60, 450 and 850/im. These 
fluxes were then subjected to 'observational uncertainties' 
by drawing a new flux from the Gaussian distribution with 
the original flux as the mean and la uncertainties derived 
from the observations as 7, 25 and 14 percent at 60, 450 and 
850/im respectively. Next, only those new fluxes which were 
found to be > 3 times the error were selected, to mimic de- 
tection limits. The required flux ratios S60/S450, Seo/Ssso 
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Figure 6. Colour-colour plot for the objects in the IRAS sample with 450^m fluxes (stars and circles), other objects from the literature 
including Mrk 231 (triangles) and some of the optically selected sample (diamonds). The conservative uncertainty in the ratios is shown 
based on flux error of 10% for 60/^m, 20% for 850//m and 30% for 450//m. The best fit line is shown, which has a slope of unity. Data for 
the stars and open circles is from this paper, triangles from Rigopoulou et al. (1996); Alton et al. (1998); Prayer et al. (1999), and the 
diamonds are from Dunne (2000) 



and S'45o/>S'85o were then calculated from these final fluxes. 
In order to be able to compare such things as the linear fit- 
ting parameters to those derived from the actual data, we 
needed to use the same number of data points (37). There- 
fore fluxes were created in sets of 37, and each set was fitted 
in exactly tfic same way as the real data (simple least squares 
analysis) and the standard deviation on the ratio 5450 /•S'gso 
calculated. 10,000 versions of Fig. 6 were created for each 
SED model, and used to produce an estimate of the mean 
slope and intercept of the best fitting line and also CT450 /gso 
(the dispersion on the ratio S'45o/'S'85o)- Additionally, one 
point firom each set of 37 was taken and placed into a sep- 
arate group. These 10,000 flux ratios were then compared 
with the distribution of the real flux ratios via Kolmogorov- 
Smiruoff two-sample tests. These points arc plotted against 
the real data for the different models in Fig. 7. We tried dif- 
ferent models to see what combination of parameters best 
reproduced the statistical properties of the data, in partic- 
ular: 

(i) the mean values and distributions of the three flux 
ratios. 

(ii) the best fitting slope and intercept. 

(iii) the dispersion of the flux ratio 5450 / Ss5o ■ 

The models represent a range of 'realistic' isothermal and 
two-component SEDs of the form fitted to the fiuxes ear- 
lier. One additional model was included to account for the 
possible contribution of the ^^CO(3-2) line to the 850/ttm fil- 
ter. This line is excited in warm {Tu > 33 K), dense gas and 
may comprise a non-negligible fraction of the bolometer flux. 
Observations of this lino arc still relatively rare as they are 
hampered by the same difficult observing conditions which 
have held continuum observers in the dark for so long. The 
only accurate way to ascertain the contribution of the line to 



the 850/im flux is to map the same galaxies in CO (3-2) and 
compare the total flux in the line to that in the continuum. 
Failing this, rough estimates of the CO (3-2) contribution 
may be made in a variety of ways (see Dunne 2000 for de- 
tails, also Papadopoulos & Allen 2000; Israel et al. 1999; 
Bianchi et al. 1998; Mauersberger et al. 1999), all consistent 
with approximately 4-30 percent of the total continuum flux 
being due to CO(3-2) line emission. This large range reflects 
the strong dependence of CO(3-2) on the physical condi- 
tions in the gas. CO (3-2) observations of 50 galaxies from 
the SLUGS IRAS sample are in progress which should al- 
low us to make a better estimate of the contamination in 
the future (Seaquist, private communication). The efi^ects of 
possible line contamination are not very large (as the sus- 
pected line contribution is comparable to the uncertainty 
in the submm fluxes) but would lead to us overestimating 
the dust masses by up to 30 percent and underestimating 
/? by ~ 0.1 in both the isothermal and two-component tem- 
perature cases. To mimick a CO (3-2) contribution Model 7 
simply adds 10 percent to the predicted 850/itm fiux from 
the SED template, which has the effect of reducing the flux 
ratios for a given /3. The tightness of the correlation in Fig. 6 
suggests that the contribution by this line adds little to the 
scatter and is therefore either relatively small or a remark- 
ably constant fraction of the continuum flux. The results of 
the comparisons of simulations and data are in Table 9. 

The two isothermal SED models (1 & 2) are not com- 
patible with the data as firstly, the ratio of ^450 / •S'sso is not 
predicted to be constant (slope not unity) , and secondly the 
means and distributions of the flux ratios are not compat- 
ible with the data and are ruled out by the KS tests at a 
high level of signiflcance. Not all two-component models can 
reproduce the observed Fig. 6 either. Models 5 & 6 which 
have lower /3 values (either /3 = 1.5 or /3 has a range of values 
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Table 8. Parameters for SED templates. 



Model Ty, (K) Tc (K) /3 Nc/N^ 



1 Uni 25 - 65 - Uni 1-2 

2 Gau n = 35.7, ct = 5.3 - Gau fj. = 1.3, a = 0.2 

3 Uni 30-55 Gau /i = 20, ct = 2.5 2 Uni 1-100 

4 Uni 30-55 18 2 Uni 1-100 

5 Uni 30-55 Uni 15-25 1.5 Uni 1-100 

6 Uni 30-55 Gau ^ = 20, ct = 2.5 Uni 1.5-2.0 Uni 1-100 
y*CO Uni 30-55 Gau ^ = 20, ct = 2.5 2 Uni 1-100 



(1) Model number. (2) Warm temperature distribution. Uni = a uniform random selection between the two limits, Gau = drawn from a 
Gaussian distribution with mean /j, and dispersion a. (3) Cold temperature distribution. No entry indicates an isothermal model. (4) 
Emissivity index. (5) Distribution of ratio of cold-to-warm dust. 

Model 2 uses the mean values and dispersions in T^j and /3 from the isothermal fits in Paper I. 

Model 7*CO [jg^g j-jjg SED parameters listed but then adds a further 10 percent to the predicted 850/im flux to mimic a contribution 

from the CO (3-2) line. 



Table 9. Results of simulations. 



(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 

Model Slope Int. ^450/850 (S60/S4B0) KS (Sso/SsBo) KS (546o/S86o) KS 



Data 
1 


1.01 ± 0.03 

0.915 
3.2a 


-0.909 ± 0.045 

-0.684 
5a 


1 6+°''2 
^•"-0.24 

2.67 
2.5e-4 


13.9 ±2.1 
64.6 


0.52 
2.4e-9 


104 ± 14 
574 


0.49 
1.9e-8 


7.90 ± 0.26 
8.13 


0.15 
0.33 


2 


0.924 
2.9a 


-0.680 
5.10- 


2.01 
0.079 


12.9 


0.1 
0.83 


89.9 


0.13 
0.51 


6.80 


0.32 
5e-4 


3 


0.993 
0.6a 


-0.910 

0.02a 


2.36 

4.1C-3 


14.6 


0.16 

0.26 


124.4 


0.09 

0.89 


8.72 


0.23 

0.03 


4 


0.987 
0.8a 


-0.881 
0.6a 


2.22 
0.0137 


15.9 


0.18 
0.18 


133.2 


0.14 
0.46 


8.34 


0.17 
0.23 


5 


0.992 
0.6a 


-0.774 
3.0a 


1.73 
0.521 


5.5 


0.5 
8.7e-9 


34.2 


0.55 
2e-10 


6.32 


0.42 
2e-6 


6 

7*co 


0.970 
1.3a 
0.992 
0.6a 


-0.808 
2.2<T 

-0.869 
0.9a 


2.14 
0.026 
2.15 
0.022 


9.0 
14.5 


0.33 
4.3e-4 
0.16 
0.26 


68.3 
113.3 


0.33 
5.7e-4 
0.13 
0.56 


7.47 
7.92 


0.20 
0.089 
0.12 
0.69 



(1) Model name, first row is the actual data. Uncertainties in the values are from the least-squares linear fitting and from the observed 
dispersions in the flux ratios. The quoted error on o"45o/g5o is a 95% eonfldence interval. The errors on the flux ratios are the errors on 
the mean value. (2) Gradient of the relationship log(S'6o/45o)i log(5'6()/85o) in Fig. 6. Determined by simple linear regression (least 
squares) on sets of 37 points. (3) Intercept of the same relationship. For the models, the value is the mean for the 10,000 simulations. 
(4) Dispersion in the ratio of 5450 /5850. (5,7,9) Overall mean ratio of fluxes (10,000 points). (6,8,10) Result of a two sample KS test on 
the 37 data points and the 10,000 simulated ratios. 

The significance of the differences between models and data are listed beneath each entry, either as no. of S.D. between means for the 
fitted slopes and intercepts, as the probability that the two samples have the same variance for a^^Q /gso (calculated from the 
F-distribution with 36 and 9999 degrees of freedom for the two samples) , or for the KS tests, as the probability that the two 
distributions are from the same parent population. 



between 1.5-2.0) are also not compatible because their flux 
ratio predictions arc too low (i.e. too much submm flux com- 
pared to that at 60^m) , however they do produce the correct 
slope. The three models which had /3 = 2 (3, 4 & 7) all pro- 
duced acceptable matches to the data (the only marginal 
point being the ratio of 5450/5850 in Model 3 which had 
a KS signiflcance of 0.03). Lowering the cold temperature 
slightly as in Model 4, solved this discrepancy. Also, Model 
7, which includes the effect of a CO(3-2) contribution can 
produce the correct flux ratio with the same cold tempera- 
ture as Model 3 (i.e. 20 K). While there is an element of Tc, 



fH degeneracy in fitting the ratio 5450/^850 (we could raise 
the cold temperature and lower /3 to produce the same ra- 
tio), this does not solve the problem with the IRAS/submm 
ratios which were all significantly different from the data in 
the lower /3 models. For this reason, we believe that Fig. 6 is 
evidence in favour of a steep and universal /3 value close to 
2. Any CO(3-2) contribution will only strengthen this con- 
clusion and further reject lower /3 values. The warm temper- 
atures used in the models were based on the values which 
came out of the two-component fits, and were kept constant 
for all models. We did experiment with altering them but 
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(a) Model 1 



(b) Model 2 



(c) Model 3 






(d) Model 4 



(e) Model 5 



(f) Model 6 




(g) Model 7 

Figure 7. Comparison of the 10,000 simulated flux ratios and the actual data points for the 37 galaxies (red stars) for models 1-7. 



for all realistic scenarios the warm temperature values did 
not have a significant efltect on the models. The same is true 
for the ratio of cold/warm dust, which was allowed to vary 
between 1-100 as per the two component fits. The tightness 
of the ratio 5450/5850 was very difficult to reproduce with 
the simulations, and indeed it was not obvious which fac- 
tors directly affected this. In all cases the models could not 
reproduce as tight a distribution as the actual data. This 
may be because we have been too conservative with our es- 
timates of the observational uncertainties, or more hkely, 
because there are intrinsic relationships between such pa- 



rameters as Tjn, Tc and Nc/N^ which were not accounted 
for in our models. 



A cold component temperature between 15-25 K and 
/3 ~ 2 is in good agreement with the measurements of 
the difltuse ISM in the Galaxy (T ~ 17 K, /3 = 2.0) by 
COBE/FIRAS (which had enough spectral resolution to do 
it properly), and with calculations of the predicted equilib- 
rium temperature for grains in the Galactic ISRF which is 
25-15 K, depending on Galactocentric distance (Cox, Kriigel 
& Mezger 1986). 
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The Verdict 

From SED fitting alone, it is often not possible to si- 
multaneously determine (3 and to show conclusively that 
there is a cold dust component, even when as for Arp220 
(Fig. 3), there are measurements at many FIR/submm wave- 
lengths. The results of the simulations of the 450/850/im 
ratio strongly favour (5 — 2, and for the 7 galaxies where a 
cold component is essential (Section 3.3), this is also the /3 
value from the fit. We will therefore assume that /3 = 2 is 
the true value and that it holds for all galaxies. Since there 
is abundant evidence in the literature for both cold compo- 
nents and /3 ~ 2 we accept it as the best physical model for 
the SED. Acceptable two-component fits using (3 = 2 were 
possible for all of the sources. 



4 DISCUSSION 

The galaxies wore re-fitted using Eqn. 5 but now assuming 
a fixed (3 = 2. The results are given in Table 10 and some 
examples of the fits are shown in Fig. 8. 



Gas and dust masses 

The dust mass has to be rc-calculatod if a colder compo- 
nent is present, and will always be higher than if only one 
temperature is assumed. The two-component dust mass is 
calculated as follows: 



Md2 



'*d(850) 



+ 



B(850/im, Tc) B(850At?n, T„) 



(8) 



where Kd(860) is assumed to have the same value as in 
Paper I (0.077 m^ kg~^). It is the temperature of the colder 
component which is critical for an accurate determination 
of the dust mass as it is this dust component which pro- 
duces most of the emission at 850/im. If Tc has been overes- 
timated then the dust mass will be underestimated, in much 
the same way as if only a single temperature had been as- 
sumed. The FIR luminosity can easily be integrated under 
the new two-component SED and this and the dust masses 
are given in Table 10, along with the revised gas-to-dust rar 
tios. On average, the dust masses increase by a factor ~ 2 
compared to the single temperature estimates from Paper 
I. In contrast, the values of Lflr do not increase by much at 
all and this is simply because the long wavelength region, 
where the cold component has modified the SED, does not 
contribute much to the integrated energy. The gas-to-dust 
ratios arc also lower than in our earlier paper by a factor 
~ 2, bringing them more into line with the Galactic value 
of 226§ (Sodroski et al. 1997). There is a slight tendency 
for the galaxies with the largest angular sizes to have higher 
gas-to-dust ratios (NGC 520, NGC 958, NGC 6181, NGC 
7541) which could be due to us missing some extended dust 
emission outside the region mapped by SCUBA. 



§ The value quoted in Sodroski et al. (1997) is 167 but we have 
had to scale tiieir kj(240) = 0.72 kg~^ to ours at 850/im. 
Using (3 = 2 the value of K(i(850) appropriate for Sodroski et 
al. is 0.057m^ kg~^. Therefore their Gj, must be multiplied by 
0.077/0.057 to be consistent with ours, giving the value 226. 




dust mass (Mg) 

Figure 9. Dust mass functions. Solid line:- 'cold dust mass' func- 
tion, dust masses arc from Table 10 and from Paper I using a sec- 
ond cold component with Tc = 20 K and (3 = 2. Dashed Unc - dust 
masses calculated from the 850^m fluxes using a single temper- 
ature as in Paper I. Schcchter parameters shown arc a = —1.85, 

a = 
dex" 



: 9.3 X 10^ Mq and (p* = 4.9 X 10" 
2.5 X 10'^ M(7) and q 



-1.23, M* 
(dashed). 



Mpc"-^ dcx-1 (solid); 
= 1.6 X 10-3 Mpc-3 



i 



'^/mMM(/mMM('/M 



cold temperature (K) 

Figure 10. Distribution of cold temperatures using (3 = 2. 

For the galaxies without useful 450/um data, one can 
obtain reasonable estimates of dust mass by fitting a two- 
component model with Tc and /3 fixed at 20 K and 2. This 
is actually the procedure used in Paper I to investigate the 
possibility of cold dust. The dust mass function for the whole 
sample using the dust masses derived in this paper and those 
'cold dust masses' (Md2) given in Table 4 of Paper I is shown 
in Fig. 9. This function could still be an underestimate if 
there are large numbers of cold, dusty galaxies which are 
not present in the IRAS sample. 



4.1 Dust heating 

The new distribution of Tc is shown in Fig. 10 and is now 
much tighter than when (3 was not fixed; it is also more 
Gaussian with a single peak at 21 K. The mean is 20.9 K 
with a = 2.5, and both the mean value and the small spread 
agrees with the range of values for Tc found in the literature 
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(a) ugc903: T(45,21), n=91 



(b) ngc958: T(44,20), n=186 



(c) ugc2369: T(36,21), n=6 





(d) ugc2403:T(50,22), n=100 



(e) irl211: T(49,24), n=27 



(f) ngcl667: T(28,17), n=3 






(g) ugc2982; T(31,17), n=6 



(h) irl525; T(45,19), n=15 



(i) ugc5376: T(44,21), n=87 






(j) ngc2623: r(50,27), n=18 



(k) arp220: T(48,18), n=42 



(1) ngc7592: T(39,19), n=20 



Figure 8. Two-component SEDs assuming fS = 2 for IRAS sample. The solid lines represent the composite 2-component SED while the 
dot-dash lines show the warm and cold components. The caption to each figure gives the two temperatures plus the ratio of cold:warm 
dust (n). 
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Table 10. Temperatures and dust masses for two-component SEDs assuming (3 = 2. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Name 






Nc 




Rm 


Lfir 


Rl 


'3d(HI) 


G'd(H2) 


Gd 




(K) 


(K) 




(log M0) 




(log Lq) 










UGC 903 


45 


21 


91 


7.37 


1.95 


10.34 


1.17 


141 






NGC 520 


50 


24 


61 


7.46 


1.74 


10.78 


1.12 


204 


234 


438 


NGC 958 


44 


20 


186 


8.28 


1.86 


11.04 


1.17 


204 


85 


289 


UGC 2369 


36 


21 


6 


8.08 


2.04 


11.34 


1.02 




269 




UGC 2403 


50 


22 


100 


7.57 


2.04 


10.71 


1.17 








UGC 2982 


31 


17 


6 


8.11 


2.34 


10.99 


1.00 


91 


132 


223 


NGC 1614 


38 


20 


3 


7.78 


2.14 


11.32 


1.00 


58 


174 


232 


NGC 1667 


28 


17 


3 


7.91 


1.82 


10.81 


1.00 


69 


117 


186 


NGC 2623 


50 


27 


18 


7.59 


1.58 


11.36 


1.07 


46 


240 


286 


NGC 2856 


45 


22 


58 


7.08 


1.91 


10.23 


1.12 


117 






NGC 2990 


42 


21 


58 


7.33 


1.91 


10.35 


1.12 


229 






UGC 5376 


44 


21 


87 


7.11 


1.91 


10.05 


1.12 


200 






NGC 3110 


42 


23 


44 


7.93 


1.58 


11.10 


1.10 




257 




IR 1017+08 


40 


19 


8 


8.20 


2.95 


11.52 


1.10 




151 




IR 1056+24 


40 


26 


6 


8.14 


1.51 


11.79 


1.02 




166 




Arp 148 


37 


21 


12 


8.28 


2.09 


11.41 


1.05 




85 




MCG+00-29-023 


36 


20 


13 


7.99 


2.14 


11.04 


1.05 








IR 1211+03 


49 


24 


27 


8.55 


2.29 


12.05 


1.12 




123 




NGC 4418 


58 


21 


62 


7.40 


3.47 


10.77 


1.15 


22 


52 


74 


UGC 8387 


33 


24 


2 


7.92 


1.29 


11.39 


1.00 




148 




Zw 247.020 


39 


25 


6 


7.50 


1.51 


11.04 


1.02 








Zw 049.057 


44 


23 


28 


7.74 


1.82 


11.05 


1.07 




65 




1 Zw 107 


39 


20 


7 


8.23 


2.63 


11.56 


1.02 




93 




Lti 1525+36 


45 


19 


15 


8.29 


3.98 


11.69 


1.00 








Arp 220 


48 


18 


42 


8.80 


3.24 


11.95 


1.02 


48 


41 


89 


NGC 5962 


31 


18 


11 


7.49 


2.19 


10.26 


1.17 


98 


120 


218 


NGC 6052 


37 


20 


16 


7.65 


2.14 


10.72 


1.07 


170 


85 


255 


NGC 6181 


31 


19 


7 


7.47 


1.86 


10.39 


1.05 


191 


195 


386 


NGC 7541 


31 


17 


6 


7.90 


2.40 


10.77 


1.02 


195 


112 


307 


NGC 7592 


39 


19 


20 


8.13 


2.75 


11.14 


1.05 


87 


129 


216 


NGC 7679 


37 


19 


14 


7.75 


2.51 


10.81 


1.05 


112 






NGC 7714 


42 


21 


13 


7.05 


2.34 


10.39 


1.05 


151 


191 


342 
















mean 


128 


142 


253 



(1) Name. (2) Warm temperature using /3 = 2. (3) Cold temperature using /3 = 2. (4) Ratio of cold to warm dust. (5) Dust mass 
calculated using the parameters in columns (2-4). (6) Ratio of two-component mass to single temperature mass from Dunne et al. 
(2000). (7) FIR luminosity (40-1000/im) integrated under the 2 component SED. (8) Ratio of this to that calculated assuming a 
single temperature SED from Dunne et al. (2000). (9) Hi Gas-to-dust ratio. (10) H2 gas-to-dust ratio. (11) Total gas-to-dust ratio (sum 
of Hi and H2). Gas masses are those used in Dunne et al. (2000) 

(15 - 25 K:- Haas ct al. 1998,2000; Frayer et al. 1999; Braine 
ct al. 1997; Dumke ct al. 1997; Calzetti et al. 2000; Guelin 
et al. 1995; Sievers et al. 1994; Neininger et al. 1996; Alton 
et al. 1998). It is also interesting to note that alternative ap- 
proaches to the SED inversion problem presented by Perez 
Garcia et al. (1998) and Hobson et al. (1993) also require two 
or more temperature components and produce similar tem- 
peratures. When observations are made with enough spatial 
resolution, the SEDs of the central and outer regions of a 
galaxy can be fitted separately. This often shows that the 
central regions contain most of the warm dust at T > 30 
K while the outer disk regions are dominated by the colder 
(15 — 20 K) dust, the total SED being the sum of the in- 
ner and outer (Papadopoulos & Seaquist 1999; Haas et al. 
1998; Dumke et al 1997; Braine et al. 1997; Trewhella et 
al. 2000). This suggests that Nc/N^ (or the prominence 
of the cold component) is really telling us about the rela- 



tive importance of the disk/bulge dust components in these 
galaxies (similar to the cirrus/starburst components in the 
models of Rowan- Robinson & Crawford 1989). The more ac- 
tive starburst galaxies in the sample display conditions like 
the centres of normal galaxies over a much larger region, so 
that the central warmer component dominates the SED. The 
fact that the warm dust always out-shines the cooler dust 
(per unit mass) is what makes it so difficult to disentangle 
the cold component from the SED in cases where there is a 
relatively large amount of warm dust. 

The equilibrium temperature of a dust grain immersed 
in a radiation field can be expressed as a function of the 
interstellar radiation field (ISRF) as (van der Hulst 1946; 
Disney et al. 1989) 

Teq oc ISRF'/^ 

This would lead one to expect variations in Tc if some galax- 
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Figure 11. 60^m luminosity correlated with dust temperature 
(from isothermal models) in Paper I. 



ies (e.g. the more active staxbursting ones) have higher IS- 
RFs - leading to sUghtly warmer grain temperatures in their 
diffuse, cold component. The range of Tc found in these 
galaxies (17-27 K) corresponds to a factor of ~ 10 in ISRF 
intensity, although it may have been expected that the more 
active star forming galaxies in the sample would have ISRFs 
many orders of magnitude greater than that in the Milky 
Way. However, these active environments arc also localised 
and dusty, therefore a lot of the extra short wavelength stel- 
lar radiation (UV) may be absorbed near to the source, pro- 
ducing enhanced FIR radiation which is in fact the warm 
component. It is generally believed that older stars con- 
tribute significantly to the heating of the diffuse 'cirrus' dust 
component as well as some radiation from OB stars which 
leaks from the star forming molecular clouds, although the 
relative importance of OB stars as a source of heating is 
still not well determined (Bothun, Lonsdale & Rice 1989; 
Cox, Kriigcl & Mezger 1986; Waltorbos & Grccnawalt 1996; 
Boulanger & Perault 1988). If OB radiation leaked from sites 
of star formation is not an important heating source for the 
colder component then the distribution of older stars and 
diffuse dust must be responsible for any changes in Tc. We 
will leave an investigation of the details of the dust heating 
to a future paper. 

In Paper I we presented a relationship between dust 
temperature (as fitted by an isothermal model) and 60/im 
luminosity (Fig. 11). Our explanation for this was that the 
position of the 60/im flux on the Wien side of the Planck 
function makes it very sensitive to dust temperature. For a 
given mass of dust, hotter galaxies would be easier to detect 
in a 60/im fiux limited sample and therefore a dependence of 
60/im luminosity on dust temperature is observed. We can 
now investigate this further using the extra information we 
have about how the dust is split into warm and cool compo- 
nents. The correlation coefficients and fits to the following 
relationships arc detailed in Table 11. We will now use the 
whole data set of 104 galaxies. Those galaxies without useful 
450/^m data were fitted with a two-component model using 
fixed 13 = 2 and Tc = 20 K (as described in Paper I). Fig. 12a 



shows that there is no correlation between the fitted warm 
temperature and 60/^m luminosity. The quantity which is 
important now is the proportion of energy in the warm com- 
ponent, which can be described roughly as (A^c/A'™)"^ x T^. 
This has a quite a strong relationship with igo as displayed 
in Fig. 12b, although there is a lot of scatter. This relation- 
ship is not linear, having a slope less than unity. It must 
be remembered here that luminosity is a distant dependent 
quantity but the relative amount of 'warm energy' is not as 
it has not been normalised in any way to the galaxy size. Is 
there a strong connection between galaxy size (or dust con- 
tent) and Leo? In Paper I we argued that the link between 
dust temperature and Lno was not a function of galaxy size 
(i.e. dust temperature and dust mass were unrelated, which 
is still the case here). However, dust mass and Leo are clearly 
related as shown in Fig. 13a although the relationship is not 
linear either (slope is less than unity) and therefore 60//m lu- 
minosity would seem to depend on both the quantity of dust 
in the galaxy and the relative amount of energy vn the warm 
component. If the 'warm energy' quantity from Fig. 12b is 
multiplied by Md2 to produce a distant dependent variable 
similar to Leo we find a remarkably tight correlation be- 
tween this product, the 'warm luminosity', and the 60/im 
luminosity (Fig. 13b) and the slope of this correlation is 
unity. Thus the typical FIR luminous objects detected by 
IRAS do not necessarily have to have vastly more dust than 
lower luminosity ones, but they do have to have a greater 
fraction of it heated to warmer temperatures. We can re- 
peat this comparison using the 'cold energy' {Nc/Nyj) x 
and 'cold luminosity' (Nc/Nyj) x x Md2. Neither of these 
quantities correlate with Leo as shown in Fig. 14. This is 
an interesting point because if there are galaxies which pos- 
sess very dominant cold components, even though they may 
have large quantities of dust, they will not necessarily have 
large 60/Ltm luminosities. This means they may not be part of 
60/im flux limited samples and hence represent a population 
so far absent from our luminosity and dust mass functions. 
This issue must await the completion of data analysis from 
our optically selected sample, which should not be biased 
against galax:ies with large fractions of colder rather than 
warmer dust. 



5 CONCLUSIONS 

We have used our own SCUBA 450/im data and fluxes 
from the literature to produce better constraints on the 
FIR/submm SEDs of the SLUGS IRAS sample. It was often 
not possible to find unique decompositions of the SEDs using 
a series of modified Planck functions, either cold tempera- 
tures and steep /3 or higher temperatures and lower /3 being 
equally acceptable. However, by the fact that the 450/850/im 
flux ratio is remarkably constant from extreme objects like 
Arp220 to lower luminosity galaxies we concluded the dust 
emissivity index /3 could be constrained to be ~ 2. Using this 
value for l3, a colder component is required in most of the 
sample, but with a very large variation in its contribution. 
The cold component has a similar temperature in all of the 
galaxies (Tav = 21 K), in agreement with those determined 
by many other authors in the literature. The dust masses 
derived using the new temperatures are higher the previ- 
ous estimates (which used a single-temperature SED) by a 
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(a) Warm temperature vs. Leo- 



(b) Warm 'energy' vs. Leo- 



Figure 12. a) Using tlie warm component temperature Tu,, tliere is no longer a relationship with Leo- b) Creating a warm 'energy' by 
combining Tyj with the ratio of cold/warm grains, there is now a correlation with Leo- 
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(a) Dust mass vs. Leo- 



(b) Warm 'luminosity' vs. Leo. 



Figure 13. a) Dust mass (using the two-component parameters) does correlate with Leo, but the slope is less than unity meaning there 
must be some other contribution to Lgo than mere size, b) Combining dust mass with the warm 'energy' to produce a distant dependent 
quantity (termed warm 'luminosity') gives a remarkably tight relationship with Leo which does have a slope of unity. 



factor ~ 2, and this brings the gas-to-dust ratios into agree- 
ment with that of the Milky Way, and other evolved spiral 
galaxies. The 60/im luminosity was found to depend on both 
dust mass and the relative amount of energy in the warm 
component, but it has no relationship with the cold compo- 
nent. Galaxies where the cold component is very dominant 
may not be well represented in GOfxm flux limited samples. 
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Table 11. Parameters for fits and correlations 



y 


X 




Prob 


a 


linear fit: y 
m 


= mx + c 
c 


Td 


log Leo 


0.44 


2.5e - 6 


4.7 


4.57 ±0.82 


-71.9 ± 19.4 


Tw 


log Leo 


0.08 


0.42 


0.8 






log [N^/N, X T6] 


log Leo 


0.62 


2.1e-12 


6.3 


0.735 ±0.068 


-9.04 ± 1.25 


log Md2 


log Leo 


0.83 


2.7e - 27 


8.4 


0.727 ± 0.036 


-9.31 ± 0.70 


log [Md2 X (Af„/Ar,)T6] 


log Leo 


0.98 


0.0 


9.9 


1.018 ±0.014 


-7.95 ± 0.29 


log [N^/N^ X T6] 


log Leo 


0.16 


0.40 


0.8 






log [Md2 X {Nc/N^)T^] 


log Leo 


0.35 


0.053 


1.9 







(1) & (2) Variables (3) Spearman rank correlation coefficient. (4) Probability that x and y are unrelated. (5) Significance of correlation 
in S.D. (6) Least-squares linear fit parameters. 
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